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Abstract. It is shown that gravitinos with mass m 3 / 2 ~ 0.1 — 1 MeV may provide suitable cold dark mat- 
\ ter candidates in scenarios of gauge mediated supersymmetry breaking (GMSB) under SO(10) grand unifi- 
■ cation coupled to supergravity, which accommodate a messenger sector of mass scale Mx ~ 10 6 GeV. This 
is due to the combined effects of renormalizable loop-suppressed operators and generic non-renormalizable 
ones governing the dilution of a pre-existing equilibrium gravitino abundance via messenger decay. The 
above range of gravitino and messenger masses can be accommodated in indirect GMSB scenarios. The 
gravitino abundance does not depend on the post-inflationary reheat temperature and it is shown that 
' leptogenesis can generate successfully the baryon asymmetry. 
Oh! 

> 

1 Introduction 

b : 

In particle physics models in which the dynamical breaking of supersymmetry is transmitted to the vis- 
ible sector through renormalizable gauge interactions, the so-called GMSB models [TJ [51 [3] , the lightest 
supersymmetric particle is the gravitino. As such it can in principle make up the cold dark matter com- 
ponent of the Universe. Previous studies have demonstrated that the messenger sector that is responsible 
for the transmission of SUSY breaking in GMSB models, may play a key role in the determination of 
the relic abundance of gravitinos [U El [7] . In particular, the lightest messenger particle with mass 
Mx > 10 5 GeV is assumed to be stable in GMSB models proposed so far. However messenger number 
must be violated in some way, otherwise this lightest messenger would overdose the Universe [5] (except 
for the particular case Mx ~ 10 4 GeV [9], [10]). On general grounds one expects this messenger to visible 
sector coupling to be suppressed in order to preserve the successful phenomenology of GMSB models, in 
particular with respect to the natural suppression of flavor changing neutral currents and correct elec- 
troweak breaking [T^]. If the decay width of the lightest messenger is sufficiently suppressed, the decay 
may be a significant source of entropy which reduces the gravitino abundance. 

Following this line of thought, we show in the present Letter that the gravitino can indeed provide 
a natural cold dark matter candidate in the "simplest" GMSB models under SO(10) grand unification 
coupled to supergravity. More specifically, we find that for a gravitino mass m 3 / 2 ~ lMeV and Mx ~ 
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10 6 GeV, generic non-renormalizable messenger-matter couplings allow to obtain the right relic abundance 
for dark matter independently of the post-inflationary reheat temperature. Leptogenesis can thus operate 
successfully at high temperatures and produce the observed baryon asymmetry. Such gravitino and 
messenger masses are predicted in indirect GMSB models such as proposed by Dine and collaborators [5] . 
In these models, dynamical breaking occurs at scale Adsb in a secluded sector and results in Mx ~ kAdsb, 
with k < 1 denoting a combination of model dependent coupling constants, and 7/13/2 ~ A DSB /mpi. Hence 
for Adsb ~ 10 6 — 10 7 GeV, one obtains Mx < 10 6 GeV and m 3 / 2 ^< 100 keV. For the same values, one 
also finds a fermion-boson squared mass splittings in the messenger sector Fx ~ M x which results in 
correct sfermion, squark and gaugino masses ~ m ~ (a/4n)Fx /Mx ~ ITeV [21 151 IT5], 

2 Gravitino abundance 

We assume that the post-inflationary Universe reheats to high temperature Trh > 10 10 GeV as gener- 
ically occurs in scenarios where the inflaton couples through renormalizable interactions to the visible 
sector. Then sparticles as well as messengers are initially brought to thermal equilibrium. Goldstinos, 
or cquivalently helicity ±1/2 gravitinos, are also brought into thermal equilibrium as the temperature 
exceeds the threshold T cq - 10 5 GeV (m 3/2 /100keV) 2 (M 3 /10 3 GeV)" 2 (M 3 gluino mass) at which pro- 
duction of goldstinos by sparticle scattering occurs faster than a Hubble time [5]. Note that gravitino 
production by messenger scatterings may further reduce this temperature |14j . Neglecting the abundance 
of helicity ±3/2 gravitinos, the number to entropy density ratio of gravitinos after reheating is thus well 
approximated by the thermal equilibrium value Y^ 2 = n 3 / 2 /s ~ 1.80 • 10" 3 (g*/230) _1 . 

If messengers sit in a pair of 16 + 16 spinor representations of 50(10), the lightest messenger is a 
linear combination of the ^-like [i.e. SU(3) x SU(2) x U(l) singlet boson] components of the 16 and 
of the 16 conjugate, as was first noted in ref. [9]. Since the singlet nature of the lightest messenger 
is a key issue for our study, it is important to qualify its genericity. Indeed, the mass degeneracy in 
these representations is lifted by various contributions to the non-singlet states: tree-level electroweak 
D— term corrections, renormalization group running effects in Mx and Fx, as well as genuine electroweak 
loop corrections. In the literature [9l [11] the focus was mainly on mass splitting in SU(2) doublets. 
Furthermore, the leading effects from the running of Mx and Fx can be milder in the GMSB models 
under consideration where these parameters are generated dynamically at scales much below Mqut- We 
have thus recomputed the loop corrections to the (squared) masses of the SU(3) x SU(2) x U(l) charged 
scalar messengers, both in the supersymmetric (Fx = 0) and susy breaking (Fx 0) configurations, 
in a minimal subtraction scheme, and taking the renormalization scale \i ~ Mx so as to resum the 
leading log effects. These corrections are found to be positive, both in the Fx — limit and in most 
of the parameter space when Fx 0, increasing the masses by ~ 10 — 20%. For instance, in the susy 
and Mx ^> mz limit, the correction to the (£>£,, Zx,)-like messenger masses is 2 sin 2g V^F while a 
somewhat larger correction obtains for the Z^-like and squark-like messengers. These corrections should 
be added to the renormalization group effects on Mx, Fx when running from the GUT scale down to the 
messenger scale. When applicable, the latter running largely dominates and increases further the masses 
of the SU(3) x SU(2) x U(l) charged messengers over that of the singlet by 20 — 60%. In the sequel we 
denote by X this lightest singlet messenger. 

The relic abundance of this lightest messenger is determined as usual by annihilation freeze-out from 
thermal equilibrium (provided the decay width to the visible sector is sufficiently suppressed, as will be 
the case here). Since X is a singlet under the standard model gauge group, it can annihilate either i) at 
the tree level into visible sector particles and only through the exchange of a GUT mass 5*0(10) gauge 
boson, ii) at the tree-level into a pair of goldstinos, Hi) at the one-loop level where the heavy as well as 
the light particles of the visible (and spurion) renormalizable sector of the model contribute in the loops. 

The contribution from i) is highly suppressed by a factor (Mx /AfcuT) 4 and can be safely neglected. 
Annihilation into goldstinos ii) is induced by supergravitational interactions (see also [8]). It receives 
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contributions from various sectors including the purely gravitational s— channel exchange of gravitons, 
4— leg contact interactions between X and the gravitino, t— channel exchange of the fermionic partner of 
X , as well as s— channel exchange of the scalar partners of the fermions which make up the goldstino after 
supersymmetry breaking. The coupling of the latter scalars to the gravitinos is fixed by the super-Higgs 
mechanism, 15, 16 , while their coupling to X has also a renormalizable non- hidden contribution through 
the coupling of the spurion to X in the GMSB superpotential. In contrast to the case of fermions or 
gauge bosons where the scale of unitarity violation is reduced from the Planck scale to much lower scales 
by the smallness of the gravitino mass [15j [16j E] , the cross section for XX — > GG has a mild high 
energy behavior, barring non-minimal Kahler contributions which are Planck suppressed. Furthermore, 
we consider the configurations where the spurion is heavier than the lightest messenger scalar, which is 
indeed the case for parts of the parameter spaced In the opposite case a quick annihilation of the lightest 
messenger to spurion pairs would occur through tree-level diagrams [accompanied by one-loop induced 
spurion decays mainly into gluons] leading typically to a too low X relic abundance. In the limit of very 
heavy spurion we find for the X particles annihilating at rest (^x*^gg") — {^-/^)^x-^xl wnere 
F = \/3m-3/2 m Pi is the goldstino decay constant. 

The contribution from Hi) to the annihilation into standard model particles is fully controlled by 
the renormalizable sector, but can be comparable to the latter in the range of parameters we consider. 
Though a detailed description of these effects is beyond the scope of the present paper we point out here 
some qualitative features in the case of the leading annihilation into two gluons. In this case there is no 
tree-level annihilation of the type i) and all contributions proceed via one-loop diagrams. They originate 
either from the spurion/messenger sector or from the SO(10) gauge sector. In the latter case, only heavy 
particles will contribute to the loops since the lightest messenger is an MSSM singlet. The GUT scale 
contributions will thus decouple leading to negligible effects, unless the pattern of 5*0(10) breaking down 
to the MSSM contains relatively light intermediate gauge sectors [18]. We will assume here that such 
patterns do not occur and neglect the effects of this sector altogether in the sequel [a detailed study of 
possible effects will be given elsewhere] . 

The main contributions to Hi) will come from the spurion/messenger sector with s-channel exchange of 
scalar spurion and one-loop vertices, and box diagrams with virtual exchange of messenger and spurion 
fermions. The renormalizability of this sector, together with the fact that both the spurion and the 
lightest messenger are SU(3) x SU (2) l x U(1)y singlets, garantees the cancellation of all infinities as well 
as the renormalization scale dependence. The resulting cross-section scales as (ciioopv) ~ / (a^ / '4tt) 2 A 4 / s , 
where A is the spurion-messenger-messenger coupling in the GMSB superpotential, the strong coupling 
constant and / is a prefactor which embodies the model dependence: a very small / occurs in the 
decoupling regime Ms 3> Mx (/ — ► when M$/Mx — > oo), but a large value is also possible close 
to resonance when Mx — Ms/2. However, more generally / ~ 0(1) for Ms ~ Mx away from the 
resonance. As we will stress later on, in the regime where (ciioop^) is very small the lightest messenger 
freezes-out while relativistic and the relic density becomes independent of the details of the cross-section. 

In the region of parameter space we are interested in, namely 1713/2 < lOMeV and Mx < 10 7 GeV, 
the annihilation into goldstinos can be neglected with respect to the one-loop annihilation cross-section 
provided / > 6 • 1CT 10 (Af x /10 6 GeV) 6 (™ 3/2 /l MeV)- 4 A^ 4 , where we have defined A .i = (A/0.1). Note 
that in this region of parameter space, the estimate of annihilation cross-section into a pair of goldstinos 
remains perturbativily reliable. This is contrary to most scenarios of gravitino dark matter in SU(5i) 
grand unification where the solution lies in a region where the cross-section seems to violate unitarity (i.e. 
(TV > 8w/M x as v — ► 0), so that the conclusions depend strongly on what one assumes for the behavior 
of this cross-section, see [8], and furthermore multi-goldstino production can then become important to 
consider. 

For instance, in the simplest model with an 5(7(6) X 1/(1) X U(l) m gauge group in the dynamical breaking sector [2], 
such configurations are achieved for V3 < ft/A < 2.2 and Ai/A <C 1, where A, ft and Ai denote respectively the spurion 
self-coupling, its coupling to the messenger fields, and its coupling to the U(l) m charged fields in the superpotential. Wider 
ranges are also found for larger values of Ai/A(< 1). 
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The relic abundance of the lightest messenger can thus be written as: 



10 6 GeVy \a 3 J 01 / 



with x f = log Q f /y/log(Qf) and Q f ~ 1.5 • 10 6 /(a 3 /4vr) 2 A^ ^(Mx/10 6 GeV)" 1 . The parameter a;/ 
denotes the ratio of the messenger mass to the temperature at which freeze-out of annihilations occurs. 
Hence 

Y X ~ 5.0 • lO- 8 /" 1 ^ (y^v) [l + 0-26 log(/) - 0.26 log (J^) + 1.05 log (A .i)] . (2) 

When 10~ 2 < A a *(M x /10 6 GeV)// < 1, Y x ^ 5 • 10- 5 (M x /10 6 GeV) a8 /-°- 8 A a f 2 provides a good 
approximation. For other values of Mx and /, one needs to keep the logarithmic correction. 

Eventually, X will decay to visible sector particles, otherwise it would overdose the Universe. As 
mentioned earlier, we assume that X decay can occur through non-renormalizable operators. The main 
motivation for this is to preserve the phenomenological successes of GMSB models, notably with respect 
to flavor changing neutral currents and electroweak symmetry breaking. It is certainly possible to intro- 
duce renormalizable operators that violate messenger number and yet do not spoil the features of GMSB 
models, however such operators have to be suppressed by unnaturally small numerical prefactors [12] . 
particularly for small messenger masses. Non-renormalizable higher dimension operators are expected 
to occur at the Planck scale, which respect the gauge symmetry of 50(10) as well as R— symmetry and 
possibly other global symmetries. Depending on the charge assignments [which should forbid dangerous 
couplings, e.g. between messengers and Higgses with GUT vev's], some of these operators will violate mes- 
senger number by one unit. For these operators, the typical decay width of X is Tx — (l/167r)/'M^ r /mp 1 , 
with /' a numerical factor of order unity. Decay occurs when the Hubble rate H ~ Tx > or equivalcntly, 
at background temperature 

/ M \ 3/2 / > \ -1/4 

The decay width of the lightest messenger is so suppressed and its relic abundance is so large (due to its 
suppressed annihilation cross-section) that X actually comes to dominate the total energy density budget 
before decaying. This happens since X quanta become non-relativistic at temperatures < Mx and there- 
fore their energy density redshifts less fast than that of radiation. This era of non-relativistic matter dom- 
ination starts at background temperature T dom ~ (A/3)M X Y X ^ 67 GeV /-°- 8 A a 5' 2 (M x /10 6 GeV) 18 , 
which indeed exceeds T > . The decay of this lightest messenger thus results in a significant amount of 
entropy generation, by reheating the Universe to temperature T > . The amount of entropy produced can 
be written in a rather simple way as [§]: 

10 2 /-°- 8 A-f/'-V 2 ( I? ^ v ) ' 3 , (4) 

and all species are diluted with respect to the background entropy density by the factor AS*. 
In particular, the goldstino abundance after lightest messenger decay reads: 



This gives a present-day relic abundance 

3.8^3.2 f/l/2 ( 

10 6 GcV/ Vo.lMeV 



^ * 0.067/" *Kf* I -°" (JgBJ) . (0) 
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Hence one finds the correct relic abundance for dark matter for m 3 / 2 ~ 0.1 — 1 MeV, Mx ~ 10 6 GeV and 
^01/ ~ /' ~ ^(1)' K the one-loop annihilation cross-section prefactor Aq x f takes values significantly 
larger than unity, the region in which one finds suitable gravitino dark matter shifts down to smaller m 3 / 2 , 
and remains at the same values of Mx- Indeed one must recall that the lower bound Mx > 10 5 GeV 
is imposed by the phenomenology of GMSB models. Note also that as Mx increases above ~ 10 7 GeV, 
depending on Aq 1/ and m 3 / 2 , annihilation into a pair of goldstinos may come to dominate the one-loop 
annihilation channel. Due to the strong dependence of the former on Mx, the relic abundance of X 
quickly decreases in this case, as Mx increases, hence entropy production becomes less effective and 
gravitinos tend to overdose the Universe. For Aq if ~ 10 3 , the region in which one finds gravitino dark 
matter lies at TO3/2 ~ 1 keV and Mx ~ 10 6 GeV. 

Conversely, if Aq x / is significantly smaller than unity, more precisely Aq 1 f < 10 _1 (Mx/10 6 GeV), 
freeze-out of lightest messenger annihilations can occur as early as in the relativistic regime, i.e. Xf < 1.5. 
In this case, the relic abundance Yx ~ 1.2 x 10 -3 does not depend anymore on the annihilation cross- 
section, and the relic gravitino abundance reads: 

n 3/2 h> * 0.028 r /2 ( f x ) f^M • (7) 

3/2 J V!0 6 GeVy VI MeV J v ' 



One may still find satisfactory solutions for gravitino dark matter, if TO3/2 > 1 MeV, Mx > 10 6 GeV 
and/or /' > 1. Note that m 3 /2 < 10 MeV is imposed by big-bang nucleosynthesis constraints on high 
energy injection if the NLSP is a bino [5]; for a stau NLSP, the constraint is generally weaker, m 3 / 2 < 
5 GeV [19]. 

It is safe to ignore the production of goldstinos during this second stage of reheating since the cor- 
responding temperature T > given by Eq. ((3]) is well below the sparticles masses. Goldstinos are also 
produced by the decay of the next-to-lightest supersymmetric particle (NLSP), which depending on 
the underlying mass spectrum, may be generally a bino or a stau. NLSPs result not only during a 
thermal freeze-out but possibly also during the messenger decays themselves. The former contribu- 
tion to goldstinos is negligible due to the entropy release that follows NLSP freeze-out; the latter is 
also generally small when compared to the previously existing gravitino abundance. The decay width 
of the NLSP is Tnlsp - (l/487r)M NLSP /( 

TO 3/2 m Pi)' so decay occurs at background temperature 

Tnlsp - 5MeV (M NLSP /100GeV) 5/2 (m 3/2 /l MeV) _ \ If one assumes that X can produce NLSPs in 
its decay, the NLSPs produced have time to annihilate before decaying to goldstinos, depending on the 
comparison of Tnlsp and T > . To be conservative, one may assume that JVnlsp NLSPs are produced per 
messenger decay, with [5]: 

N „ M * „ * fa/2 ( Mx \ 5/2 / Mnlsp \ 2 

If these can decay to goldstinos directly without annihilating, the amount of gravitinos produced is: 

v nlsp AT Y x m -7 f i ( M x \ ( Mnlsp \ 2 , Q s 
Y 3/2 -V NLSP — -2-10 f[ W Sev) {-^Gev) - ( 9 ) 

which remains small compared to Y3/2 calculated in Eq. ((5]) above provided Mx < 10 6 GeV and/or 
Mnlsp <L 100 GeV. Moreover NLSP annihilations prior to NLSP decay, reduce this estimate [20j [5] to 

yNLSP^.o in -ii./-i/af ( g NL S p^) y 1 f M x V 3/ V g> V /4 nn , 
Y ^ - L8 ' 10 f Uo- 7 GeV- 2 J [WGev) {mjEj ' {W) 

provided Tnlsp < T y , which is the case unless /' is very small. Annihilations are quite effective for stau 
NLSPs with typical (ctnlspw) ~ 10- 7 GcV" 2 (m N LSp/100GeV)- 2 . 
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The present scenario respects the various constraints from big-bang nucleosynthesis (BBN). For in- 
stance, the lightest messenger and the NLSP decay well before nucleosynthesis provided Mj > 10 5 GeV 
and TO3/2 < 10 MeV so that their decay products can thermalize and constraints from high energy injection 
do not apply. On the other hand, for gravitino masses ~ 0.1 — 1 GeV decay occurs during BBN, and for the 
right NLSP abundance and hadronic branching ratio, interesting effects on the 7 Li and 6 Li abundances 
may occur [21] . Nevertheless, such high TO3/2 are successfully accommodated in the present scenario only 
for /'< 1. The gravitinos produced in NLSP decay, although they remain highly relativistic at the time 
of big-bang nucleosynthesis, are in too small numbers to contribute significantly to the energy budget. 
One indeed evaluates the fractional contribution to the radiation energy density carried by those relativis- 
tic gravitinos at BBN to be - F 3 ^ 2 jSP M nl sp/T nl sp ~ 2 • 10 4 y 3 ^ SP (M NL s P /100 GeV)- 3 / 2 (m 3/2 /l MeV), 
hence it can be neglected. If R— parity does not hold, the NLSP can decay much earlier through other 
channels, and BBN constraints are actually replaced with constraints on the gravitino lifetime from 
the non-observation of cosmic diffuse backgrounds distortions. However, it has been shown that if 
TO 3/2 10 MeV, these constraints are satisfied [32]: the gravitino is then so long-lived that it can be 
considered as stable on our cosmological timescale. Hence the present scenario for gravitino dark matter 
remains valid even if R— parity is violated. 

Finally, note that the bulk of gravitinos produced here behave as cold dark matter from the stand- 
point of large scale structure formation, even for gravitinos masses as small as ~keV, due to the cooling of 
gravitinos during entropy release. One indeed calculates [5] that the present-day gravitino velocity vq ~ 
0.0017km/s(m 3/2 /lkeV)- 1 / - 26 Ao:? 4 / /1/6 (M x /10 6 GcV)'- - 1 for the dominant component that results 
from thermal equilibrium (which excludes the non-thermal part from NLSP decay). The corresponding 
smoothing spatial scale for structure formation can be calculated as R ~ 235 kpc (?;o/0.05kni/s) 086 [23"] , 
which confirms that these gravitinos are approximately cold. This situation is contrary to that encoun- 
tered in scenarios in which more massive dark matter gravitinos m 3 / 2 ^ 10 — 100 GeV are produced by 
non-thermal processes (in particular by the decay of the NLSP); in these models, indeed, the gravitino 
dark matter is warm in a large part of parameter space [24] . 

In conclusion, for a reasonable choice of parameters of the underlying GMSB model, one can obtain 
the right amount of gravitinos to explain the cold dark matter content of the Universe. The vast majority 
of these gravitinos have been produced in scatterings (or decays) at high temperatures and cooled and 
diluted by the lightest messenger out-of-equilibrium decay. A significant advantage of the present scenario 
is that the final abundance of gravitinos is independent of the post-inflationary reheating temperature 
Trh- In particular, the stringent constraints on Trh for light gravitinos [T41I25] are irrelevant here. This 
implies notably that scenarios of leptogenesis from right-handed (s)neutrino decay can operate at high 
temperatures and produce the observed baryon asymmetry. This possibility has been discussed in some 
detail in a related context by Fujii & Yanagida [5J. Here we omit the details for simplicity. The decay of 
each right-handed (s)neutrino yields a net lepton asymmetry: 

|£l1 ^li^^' 

with Mm the lightest RH (s)neutrino mass, m„3 ~ 0.06 eV the heaviest left-handed neutrino mass, (H®) 
a Higgs vev, and <5 e ff an effective CP— violating phase. Assuming that the RH (s)neutrino is initially in 
thermal equilibrium (if Trh ^ Mm), the net baryon asymmetry produced reads [5]: 

^«3.6-10- 3 -^| ei |a, (11) 
s AS 

where C = —8/23 in the minimal supersymmetric standard model denotes the effectiveness of L to B 
conversion, and a ~ 1 characterizes the fraction of lepton asymmetry surviving after RH (s)neutrino 
decay. One finally obtains: 
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nB ~ I a 10 -io f o.a A a.2 f /i/af M x \ 3 ( M Rl \ 
s 1 Aoi/ Uo 6 GevJ U-10 n GevJ' [U) 

for 5 e ff ~ a ~ 1, which matches the measured asymmetry ~ 8 • 10 _ . 



3 Conclusions 

We have shown that a gravitino of mass m 3 / 2 ~ 0.1 — 1 MeV provides a natural cold dark matter candidate 
in 50(10) GMSB scenarios coupled to supergravity, with messenger mass scale Mj ~ 10 6 GeV and non- 
renormalizable messenger-matter coupling. Although our conclusions are largely independent of the 
details of the underlying GMSB model, we find that our scheme can be accommodated in the so-called 
indirect "simplest" GMSB models a la Dine et al. [2J, in which the dynamical supersymmetry breaking 
in a secluded sector is fed down radiatively to the messenger sector. For a reasonable choice of the 
various model parameters, e.g. the coupling constants, one obtains the correct dark matter abundance. 
Our results, and in particular the required combination of TO3/2 and Mx, depend on the numerical 
prefactor / of the one-loop annihilation cross-section of the lightest messenger, which is a £{7(3) x 
SU(2) x U(l) singlet scalar, the spurion-messenger coupling constant A, as well as on /' the prefactor of 
the decay width of the lightest messenger via non-renormalizable messenger-matter interactions. Non- 
renormalizable interactions are preferred over their renormalizable counterparts in order to not violate 
existing limits on flavor-changing neutral currents. On the other hand, messengers have to decay to be 
cosmologically acceptable. When /A 4 is sufficiently small, our results do not depend on the magnitude 
of the messenger annihilation cross section. 

Our results do not depend on the mass spectrum in the visible sector. In effect the final relic gravitino 
abundance is simply the abundance of goldstinos at thermal equilibrium diluted by the entropy produced 
in the lightest messenger late decay. Goldstinos are initially brought in thermal equilibrium by scatterings 
and decays involving sparticles as well as messenger fields. This fact is at variance with other SUSY 
models in which the dark matter candidate is a neutralino LSP; in those models, one must find the right 
combination of the parameters that determine the visible sector mass spectrum, e.g. mi/2, too, (^0), 
tan/3 and sign(/Lt) to obtain the right dark matter abundance. 

Gravitino dark matter in this mass range 7713/2 >keV cannot be observed in dark matter search 
experiments, either direct or indirect. However there exist interesting proposals to detect evidence for 
gravitinos/goldstinos [26) in next generation colliders. Such experiments should also confirm or dispute 
the GMSB phenomenology which leads to a highly predictive mass spectrum with distinctive features. 
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